Spherical Boson Stars as Black Hole mimickers 



O 
Oh' 



F. S. Guzman 1 and J. M. Rueda-Becerril 1,2 

1 Institute) de Fisica y Matemdticas, Universidad Michoacana de San Nicolas de Hidalgo. Edificio C-3, 
Cd. Universitaria, C. P. 58040 Morelia, Michoacdn, Mexico. 
2 Facultad de Ciencias, Universidad Autonoma del Estado de Mexico, 
Instituto Literario 100, Toluca, Edo. Mex., 50000 Mexico. 
(Dated: September 8, 2010) 

We present spherically symmetric boson stars as black hole mimickers based on the power spec- 
trum of a simple accretion disk model. The free parameters of the boson star are the mass of the 
boson and the fourth order self-interaction coefficient in the scalar field potential. We show that 
even if the mass of the boson is the only free parameter it is possible to find a configuration that 
mimics the power spectrum of the disk due to a black hole of the same mass. We also show that for 
each value of the self-interaction a single boson star configuration can mimic a black hole at very 
different astrophysical scales in terms of the mass of the object and the accretion rate. In order 
to show that it is possible to distinguish one of our mimickers from a black hole we also study the 
deflection of light. 

PACS numbers: 04.70.-s 04.40.-b 05.30.Jp 
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I. INTRODUCTION 

Due to a number of observations related to high en- 
ergy events, at present time one important problem in 
relativistic astrophysics is that of the nature of black hole 
candidates (BHCs). It is usually assumed that black hole 
solutions are the only models of BHCs. This together 
with the construction of various non- vacuum solutions in 
general relativity for various types of matter and config- 
urations, the question of whether BHCs are black hole 
solutions or a different solution has been of certain in- 
terest. Among such alternative solutions are wormholes 
[l|, 0, gavastars 0, H|-[5[, brane world solutions @ and 
boson stars • The reason is that astrophysical impli- 
cations of these solutions cannot be distinguished from 
those due to black hole solutions at distances far from 
the event horizon region and thus the existence of the 
horizon itself has become the object of study of BHCs 

Various properties of black hole mimickers are to be 
studied, being in the first place the stability. For in- 
stance, stability of the whole set of wormholc solutions 
usually thought of as mimickers is not clear, instead the 
instability of wormhole solutions has been shown for ba- 
sic wormhole solutions flll - [l4| , destroying previous hope 
on the possibility that these are stable as shown in [l5[ for 
particular types of perturbations. On the other hand, the 
stability of gravastars has been explored and has been 
found that there are stability regions 
on favor of boson star solutions it can be said that the 
stability of the solutions has actually been exhausti vely 
studied, for instance using perturbative methods fl8lll9l]. 
catastrophe theory [2fj| and full non-linear numerical rel- 
ativity, both in spherical symmetry [2l|, [22j and full 
3D [23|, and the stable branches of solutions are well 
known. 

In fact, the study of boson stars has beenpushed for- 
ward up to the binary boson star collision |24j and the 



non-linear evolution of perturbed boson stars [25j . in 
both cases considering the system as a source of grav- 
itational waves. The study of gravitational wave signa- 
tures has also been studied in order to tell between a 
black hole and a gravastar Q. Instead, in the case of 
wormholes e.g. simple solutions (supported by a phan- 
tom scalar field) would not stand the perturbative anal- 
ysis and there is no hope for a study of a binary system, 
because as shown in [13[ the life time of these solutions 
is rather short and they should either collapse and form 
black holes or explode before they could merge, although 
in 1J] it was shown that charged wormholes can have a 



longer life-time if the charge parameter is adequately cho- 
sen. About a potential full non- linear study of gravastars, 
the challenges are related to the evolution of the fluid 
and numerical methods should be able to handle very 
well located distributions of matter, nevertheless some 
alternative gravastar-like solutions have been proposed 
to ameliorate the infinitesimally thin shell problem by 
introducing anisotropies in the fluid (2t|. At the end of 
the day. if black hole mimickers are to be compared from 
all the angles with black holes, mimickers should also be 
expected to do what black holes can do: they collide and 
generate gravitational radiation with a given fingerprint. 

The black hole mimicker models mentioned above are 
usually compared using the study of the accretion of ge- 
ometrically thin optically thick disk models, since they 
are based on the study of time-like geodesies on a fixed 
background space-time. In this paper we present a set of 
boson star solutions that can act as black hole mimick- 
ers based on the same accretion disk model as done for 
the other mimickers. We show that they appear to be 
independent of the astrophysical parameters related to 
the mass of the BHC and the accretion rate. With these 
results we find that the power spectrum obtained out of 
this type of accretion disk model does not suffice to dis- 
tinguish a boson star from a black hole if the boson star 
is chosen appropriately. This paper is a generalization of 
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[9j], in the following sense: 1) we show that a boson star 
mimicker for one BHC mass is also the same mimicker 
for other BHC mass by only changing the values of the 
boson mass, 2) we track the mimickers for various values 
of the self-interaction parameter and 3) we present the 
deflection of light as a method to eventually distinguish 
a black hole from a boson star. 

In order to provide potential predictions that help dis- 
tinguishing our boson stars from a black hole, we study 
the light deflection by our boson stars and by the mim- 
icked black hole, show that in the boson star case there 
is no photon sphere as expected and show the scale at 
which high resolution lensing would distinguish between 
the boson star and the black hole. 

The paper is organized as follows: in section |n] we 
present the construction of boson star solutions. In sec- 
tion III we describe the disk model and the calculation 
of the power spectrum. In section IV we present the bo- 
son stars able to mimic black holes whereas in section V 
we present a prediction to falsify our model based on the 
deflection of light by boson stars. Finally in section VI 
we draw some conclusions. 



II. SPHERICAL BOSON STARS 

Boson stars (BSs) are solutions to Einstein's field equa- 
tions minimally coupled to a complex scalar field [27l [28j . 
The scalar field is endowed with a potential of self- 
interaction that guarantees a global U(l) symmetry. This 
global phase invariance implies the existence of a con- 
served charge corresponding to the boson particle num- 
ber of equilibrium configurations. The Lagrangian den- 
sity describing a minimally coupled complex scalar field 
to general relativity reads: 
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where kq = 8?t in units where c = G = 1, <f> is the 
scalar field, the star stands for complex conjugate and 
V is the potential of the scalar field. Notice that this 
Lagrangian density is invariant under the global t/(l) 
group, and the associated conserved charge is called the 
number of particles (defined below). When the action 
is varied with respect to the metric Einstein's equations 
arise — KqT^, where the stress-energy tensor reads 



4 + d^d„cf ) *}--g^*- a cb, a +V(\cj > \ 2 )}. 

(2) 

In the case of boson stars as black hole mimickers, it 
suffices to consider the potential of the type V(|</>|) = 
m 2 \(j)\ 2 + A|^>| 4 where m is associated to the mass of the 
boson and A to the self-interaction of the boson system. 
Finally, Bianchi identities reduce to the Klein-Gordon 
equation 
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where □ </> = 

Boson stars are spherically symmetric solutions of the 
equations above when a harmonic time-dependence for 
the scalar field is assumed 4>(r,t) — o ( r ) e_Jt * ;t j where r 
is the radial coordinate and t the coordinate time. This 
condition implies that the stress-energy tensor becomes 
time- independent, which in turn implies that the geom- 
etry of the space-time is also time-independent. 

Boson star solutions are constructed considering the 
time-independent spherically symmetric line element 



ds z 



-a{r) 2 dt 2 + a{r) 2 dr 2 + r 2 dtt 2 , 
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with a(r) the lapse function and a(r) the radial metric 
function. Under these specifications, the Einstein-Klcin- 
Gordon equations become: 
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The system ([5]) is a set of coupled ordinary differential 
equations to be solved under the conditions a(0) = 1 that 
provides spatial flatness, <fio(0) finite and d r <fio(0) = as a 
condition that guarantees the regularity of the operators 
at the origin; also demand </>o(oo) = in order to ensure 
asymptotic flatness at infinity. The problem turns into 
an eigenvalue problem for w, and for a given value of 
the central field O there is a unique uj with which the 
boundary conditions are satisfied. In order to get rid of 
the constants in the equations we re-scale variables and 
constants by </>o = v/^^o, f — mr, t = u>t, a = 
and A = 2X .> . In terms of these new constants, and 
after removing the tildes from the variables the system 
of equations becomes 
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Notice that lo now turns into the central value of the 
lapse a due to the rescaling. This is the system we 
solve using finite differences with an ordinary integrator 
(adaptive step-size fourth order Runge-Kutta algorithm 
in the present case) and a shooting routine that bisects 
the value of to (the central value of a). Also, it is worth 
noticing that the radial coordinate is scaled by m and 
thus a natural scale for our calculations will be given by 
m = 1 which together with A are the two parameters 
determining the properties of each configuration. 

The solutions of ([5]) define sequences of equilibrium 
configurations like those shown in Fig. [T] Each point 
in the curves corresponds to a boson star solution. In 
each of the curves two important points for the differ- 
ent values of A are indicated: i) the critical point - 
marked with a filled circle- indicating the threshold be- 
tween the stable and unstable branches of each sequence, 
that is, configurations to the left of this point are sta- 
ble and those to the right are unstable as found for bo- 
son stars in the full non-linear regime [2p - [23j and ii) the 
point at which the binding energy Eb = M — Nm = 0, 
indicated with an inverted filled triangle, where N = 
/ j°d 3 x = J ! v /=: gg'" / [</>*<% - 4>d„(f>*}d 3 x is the num- 
ber of particles; that is, the conserved quantity due to 
the invariance under the global U(l) group of the La- 
grangian density (fTJ). For a mass estimate we use the 
Misner- Sharp mass defined through the metric functions 
by M(r) = (I — l/a 2 )r/2, being M the mass function 
evaluated at the outermost point of the numerical do- 
main; the configurations between the instability thresh- 
old and the zero binding energy point have negative 
binding energy (Eb < 0) and collapse into black holes, 
whereas those to the right of the inverted tr iang les have 
positive binding energy and disperse away [23j |, Those 
configurations to the left of the threshold of instability, 
that is, stable configurations, obviously posses negative 
binding energy, and are the expected ones to be studied 
as black hole mimickers. For reviews and more informa- 
tion on boson stars see [23, u3fj ■ 



III. DISK MODEL AND THE POWER 
SPECTRUM 



Following a procedure similar to that described in [8( 
to calculate the power spectrum from the disk model, 
we first study the geodesies. Given the line element 
ds 2 = —a(r) 2 di 2 + a(r) 2 dr 2 + r 2 dfl 2 , the equation for 
time-like geodesies followed by test particles on the equa- 
torial plane reads: 
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FIG. 1: (Top) Sequences of equilibrium configurations for dif- 
ferent values of A are shown as function of the central value 
of the scalar field </>o(0); each point in the curves corresponds 
to a solution of the eigenvalue problem and represents a bo- 
son star. The filled circles indicate the critical solution that 
separates the stable from the unstable branch. Those con- 
figurations to the left of the maxima represent stable config- 
urations. The inverted triangles indicate the point at which 
the binding energy is zero. The mass of the configurations is 
used in units of the squared Planck mass over the mass of the 
boson, which implies a scale invariance of this plot under the 
change of m. (Bottom) We show the compactness of the equi- 
librium configurations shown. We define compactness as the 
ratio 2Ngg /-R99 , where Rgg is the radius at which 99 percent 
of the total number of particles iVgg is contained. In this plot 
we also show the critical point with a filled circle and indicate 
that the critical configurations are the most compact ones 
among the stable configurations for each value of A. Those 
configurations to the right of the critical solutions are even 
more compact however unstable and are not considered here. 
In our results below we will also use the critical solutions for 
reference as the most compact stars. 



, if L 2 \ E 2 
r 2 + — 1 + — =-2-2' ( 7 ) 

where L 2 = r A ip 2 and E 2 = —a 2 i 2 are the squared angu- 
lar momentum and energy at spatial infinity, which are 
the conserved quantities of the test particle related to 
the independence on the azimuthal angle (p and t of the 



space-time respectively; an overdot indicates derivative 
with respect to the proper time of the test particle. The 
geodesies for a Schwarzschild black hole are given by (O 
with the values a 2 = a~ 2 = (I — — ). Since the equation 
for a is linear, we rescalc this function so that at infin- 
ity a(r — > 00) = l/a(r — > 00), thus at infinity a 2 = a~ 2 , 
which implies the coefficient of E 2 in ([7]) equals one for 
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both boson stars and black holes. 

The study of stable orbits of test particles requires the 

effective potential V 2 ^ = ^3-^1 + -^-^ for a given angu- 
lar momentum of the test particle. In order to consider 
only circular orbits we rewrite for r = and demand 



dV/dr = with V = [(l + L 2 /r 2 ) - E 2 /a 2 ] ft 
condition implies 
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where a 2 — raa' > for all BSs. Before proceeding to the 
calculation of the emission spectrum we need the angular 

velocity of a test particle given by il — J which is a 
regular function for all values of r in the boson star case. 

The accretion disk model is that of a geometrically 
thin, optically thick, steady accretion disk. The power 
per unit area generated by such a disk rotating around 
a central object is given by @, l3l| and used also in the 
study of disks around gravastars [5| : 
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where M is the accretion mass rate and r, is the inner 
edge of the disk. For black holes this radius is assumed 
to be at the ISCO (r = 6M) of the hole. For BSs we 
choose n — based on two considerations: i) BSs allow 
circular orbits in the whole spatial domain and thus there 
is no kincmatical restriction to choose r, = and ii) it 
has been shown that the luminosity of the disk for BSs 
considering this inner radius never reaches the value of 
Eddington luminosity in the whole spatial domain even 
for high accretion rates and thus there are no radiation 
pressure effects imposing a restriction on the inner edge 
location Q . Considering the disk starting at i\ = does 
not affect the boson star structure as long as the disk 
is made of test particles; the case where the interaction 
of the disk matter with the geometry of the space-time 
is taken into account -which is beyond the scope of this 
paper- may impose restrictions to the geometrical prop- 
erties of the disk. Moreover, considering other values of 
r,i for the BSs within the present model permits eventu- 
ally to find also another configuration that would act as a 
mimicker. Now, assuming it is possible to define a local 
temperature we use the Steffan-Boltzmann law so that 
D(r) = aT 4 , where a = 5.67 x 10~ 5 erg s" 1 cm~ 2 K~ 4 
is the Steffan-Boltzmann constant. Now, considering the 
disk emits as a black body, we use the dependence of T on 
the radial coordinate and therefore the luminosity L(y) 
of the disk and the flux Fty) can be calculated using the 
expression for the black body spectral distribution: 



Lty = And 2 F{v) 
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where d is the distance to the source, r, and 77 indicate 
the location of the inner and outer edges of the disk, 
h = 6.6256 x 10 -27 erg s is the Planck constant, 
k = 1.3805 x 10~ 16 erg K^ 1 is the Boltzmann constant 
and 1? is the disk inclination. The algorithm to construct 
the emission spectrum for such a model of accretion disk 
around a BS and a BH is as follows: 

1) Define the space-time functions a and a by choosing 
one of the equilibrium configurations in Fig. Q] and 
calculate M. 

2) Define the metric of the equivalent BH through 
a BH = v/1 - 2M/r and a BH = W 1 - 2M/r. 

3) Calculate the angular velocity, angular momentum 
and energy of a test particle for both space-times 
^bs,bh, Lbs,bh, Ebs,bh- 

4) Use such quantities to calculate the power emitted in 
both cases Dbs{t) an d Dbh(t) defined in ©■ 

5) Calculate the temperature of the disk in both cases 
T BS {r) = {Dbs^/o) 1 / 4 and T BH (r) - {D BH {r)/ a) 1 ^ . 

6) Use such temperature to integrate the luminosity (|T0|) 
Lb sty) an d L B B\ty) for several values of v. 



IV. BOSON STARS AS MIMICKERS OF 
BLACK HOLES 

We use the simple model described above and look for 
the boson star configurations that are able to produce 
approximately the same power spectrum due to the pres- 
ence of a black hole of the same mass. In order to do so, 
we first set two physical situations: 

A) a black hole of 3xl0 9 M Q with an accretion rate of 
2 x 10" 6 Af o /yr. 

B) a black hole of 10M Q with an accretion rate of 2 x 
10- 12 M & /yr. 

The angle used in the calculations is 1? = 60 degrees in 
all the cases presented, and the disk outer edge is as- 
sumed to be located at 50 M. A different inclination of 
the disk would shift the power spectrum toward higher 
luminosities for small "d and smaller luminosities for •& 
approaching 7r/2 in all the frequencies. 



A. The unexpected case of A = 0. 

In the past it was considered that BSs could be de- 
tected because the emission power spectrum from an ac- 
cretion disk showed a hardening at high frequencies by 
orders of magnitude [1], however the boson star chosen 
was the most compact in the stable branch of the A = 
case; it was shown in Q that by choosing a different bo- 
son star with a different compactness it was possible to 
mimic the spectrum. Here we show that the spectrum is 
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mimicked independently of the physical scale at the price 
of changing the value of the boson mass. 

In Fig. [3] we show the two physical situations A and 
B and mimic the black hole accretion disk spectra with 
one corresponding to a disk around an appropriate bo- 
son star configuration. For reference we also present the 
spectra calculated using the most compact configuration 
in the stable branch of the A = case (see Fig. 1) that 
corresponds to a boson star mass M = 0.633M 2 ; /to. On 
top of the black hole spectrum there is the one due to 
a boson star configuration with M = 0.473Af 2 ;/TO on 
the same stable branch of boson star configurations. In 
the case A a mass of 3 x 1O 6 M and M = 0.473M 2 j/to 
needs the mass of the boson to be m = 1.2 x 10 _27 GcV, 
whereas for the case B, a BHC mass of 10M Q needs 
to = 2.51 x 10 _22 GcV. That is, we show here that it is 
actually possible to mimic a black hole at two very differ- 
ent scales at the price of changing the mass of the boson. 
This indicates that restrictions on the value of m should 
be in turn; at the moment, since no fundamental scalar 
field particle has ever been discovered the boson mass is 
a free parameter; once there are bounds on the mass of 
a boson particle the set of BHCs that can be mimicked 
with boson stars will be restricted, and the only free pa- 
rameter will be the self-interaction. 



B. The case A / 0. 

With the introduction of the self-interaction fourth or- 
der term in the potential, it happens that a new free pa- 
rameter is involved, which makes easier to fit for instance 
the power spectrum of a disk, and is bad for the boson 
star model in terms of the number of parameters. How- 
ever, if a boson is to be observed in the laboratory, the 
mass of it would be fixed and the A parameter would be- 
come the only free parameter. Meanwhile here we present 
a few cases for A ^ and indicate the range in which the 
boson stars as black hole mimickers should be searched 
in the configurations set of Fig. 1. 

We consider again the physical situations A and B in 
order to show that the best mimicker is scale independent 
as in the case A = 0, that is, we chose our best black hole 
mimicker from Fig. 1, for a given A in one case and it 
will also work for other astrophysical cases with different 
masses of the central object and different accretion rates. 

We choose the case A = 50 and present the spectra 
as before for cases A and B. Again, there is a mimicker 
configuration with mass M = 0.9898A/ 2 ;/to. In case A 
the mass of the boson required for the boson star to be 
a mimicker is to = 1.88 x 10~ 27 Gev whereas for case B 
to = 5.25 x l(T 22 GcV. 



C. Tracking all the cases. 

The mimickers based on this model of accretion disk 
have been located in the diagram of Fig. 1 for each value 
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FIG. 2: The emission spectra for cases A (top) and B (bot- 
tom) due to the black hole and two different boson stars: 
the most compact of the stable branch for A = case 
(M = 0.633M^/m, see Fig 1.) and an adequate configura- 
tion (M — 0.473Mp;/m) that presents a power spectrum that 
mimics that of the black hole. The most compact boson star 
is used only as a reference. The configuration with M = 0.473 
is used for the two different astrophysical situations. For case 
A: m = 1.2 x KT 27 GeV and for case B: m = 2.51 x l(T 22 GeV. 
The spectrum due to the black hole and to the mimicker lie 
approximatelyon top of each other. For similar physical con- 
figurations the spectrum is comparable to those obtained in 
[a|, where also the redshift was taken into account. 



of A as those configuration which spectrum approximates 
that of the black hole with the same mass, only in terms 
of the correct boson mass. In order to show where the 
mimickers we found are located among the whole set of 
boson star solutions, we present Fig. 0] where we indi- 
cate the location of the mimickers and their apparently 
continuous distribution. 

It is natural to expect that the mimicker should change 
if the disk model is different, if the angle of inclination 
used to calculate the luminosity is different and if other 
effects are taken into account. Nevertheless, our results 
show that the boson star mimickers presented here are a 
set of measure zero among the total set of stable boson 
stars with different values of to and A out of which correct 
mimickers can be chosen. 
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FIG. 3: The emission spectra for cases A (top) and B (bot- 
tom) due to the black hole and two different boson stars: 
the most compact of the stable branch for A = 50 case 
(M = 2.257M|;/m, see Fig 1.) and an adequate configuration 
(M = 0.9898A/pj/m) that presents a power spectrum that 
mimics that of the black hole. As in the previous case, the 
same configuration with M = 0.9898 is used for the two differ- 
ent astrophysical situations. For case A: m = 1.88x 10 _27 GeV 
and for case B: m = 5.25 x 10 _22 GeV respectively. The spec- 
trum due to the black hole and to the mimicker lie approxi- 
mately on top of each other. 
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FIG. 4: The mimicker configurations are indicated with stars 
for values of A = 0, 10, 20, 30, 40, 50, 60, 70. It is expected that 
for every value of A there is a mimicker configuration that lies 
near the ones indicated for the values shown here. 



Another point has to do with the boson masses in terms 
of the mass of the BHC. We show in Fig. [5] the boson 
mass scaling with the mass of the black hole candidate for 
mimicker configurations. As an example, assuming the 
axion mass is of the order of lmeV (10 _12 GeV) there is no 
astrophysical configuration even about 1 Mq that might 
possibly act as an astrophysical black hole mimicker; in- 
stead the mass of the object would be 10~ 8 Af Q . An- 
other example is the Higgs boson, with a mass bounded 
by ftiHiggs > 115GeV would imply black hole mimickers 
with mass M ~ 10~ 22 M n . 
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FIG. 5: We show the mass of the boson in GeVs in terms 
of the mass of the black hole candidate in solar masses for 
mimicker configurations and various values of A. We extend 
the case A = case in order to show that for the axion mass 
case the object would be very small. 



LENSING, A WAY TO DISTINGUISH A 
BLACK HOLE FROM ITS MIMICKER 



In order to provide a method to distinguish a boson 
star mimicker from a black hole we show the lensing de- 
flection effects for a boson star mimicker configuration 
and the black hole mimicked. For this we choose to solve 
the full geodesies equation of null rays for the mimicker 
and the black hole 



(fix 1 ^ dx J dx 
d'y 2 J dj dj 



= 0, 



(11) 



where x l — (f , r, <p) arc respectively the time, radial, az- 
imutal angular coordinate, and 7 is the affine parame- 
ter of the geodesies. We integrate the geodesic equa- 
tions numerically using a fourth order accurate Rungc- 
Kutta integrator in terms of the affine parameter. For 
the black hole space-time we use de Schwarzschild solu- 
tion in Schwarzschild coordinates. For the boson star, 
since the metric functions are calculated as a numerical 
solution of equations ©, we solve the geodesic equation 
by sourcing the ODE integrator with second order inter- 
polated values of the Christoffcl symbols in (fTTj) . In the 
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calculations of the null ray trajectories we have verified 
the convergence of the numerical method by monitoring 
the null condition of the tangent vector of the geodesic 
at every point of the space. 

We use various values of the impact parameter and 
calculate the deflection angle in both space-times (the 
black hole and the boson star) . The results are shown in 
Fig. |5] as a set of null rays on the equatorial plane for 
the black hole and its mimickcr for the case A = 0. It 
can be seen that for small values of the impact parame- 
ter the black hole traps the null rays, whereas the boson 
star is transparent and only deflects the null rays. The 
deflection angle in terms of the closest approach to the 
center is shown in Fig. where it is shown the expected 
result that boson stars do not have a photon sphere since 
they have no horizon, that is, there is no minimum im- 
pact parameter after which the null ray trajectories are 
trapped. 

In the example of Fig. [5] the luminous source is located 
at a distance of 50M. This distance can be extended ar- 
bitrarily far away from the BHC. From Fig. 7 it can be 
learned that if the deflection angle can be measured with 
observational resolution of r ~ 15 — 20M, a boson star 
mimicker can be distinguished from a black hole because 
the curves are clearly distinct. On the other hand, be- 
cause the deflection angle due to the boson star mimicker 
starts being on top of that of the black hole for r > 20M, 
these objects cannot be distinguished from each other us- 
ing the deflection angle in this regime. 




r/M = 2 ■ 
Boson Star mimicker ■ 
Black Hole 



20 



15 



10 





x/M 



10 



15 



FIG. 6: We show a set of null geodesies for the black hole and 
its boson star mimicker. We choose the mimicker for A = 0. 
We use units for the coordinates on the equatorial plane x 
and y scaled with M (the mass of the BHC) and therefore 
this plot corresponds to both cases A and B for A = 0, where 
the scales for cases A and B are obtained when the value of 
M of the BHC is substituted. 
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Black hole - 



2 4 6 8 10 12 14 16 18 20 22 

FIG. 7: We show the deflection angle of the null rays for the 
black hole and its boson star mimicker corresponding again to 
the case A and A = presented above. In this plot r v is the 
photon sphere radius of the black hole, rs is the Schwarzschild 
radius, and the label r c indicates the closest approach radius. 
We only show the deflection angle for r > M because in spher- 
ical coordinates the geodesic equations become stiff; a detailed 
study of this region would need a patch with different coordi- 
nates which is not within the scope of this paper. 



spectrum of a simple accretion disk model around a black 
hole is compared to another around a boson star of the 
same mass. 

For each value of the self-interaction coefficient (in- 
cluding the free field case (A = 0), given astrophysical 
parameters for the mass of the black hole candidate and 
the accretion rate parameter, it was possible to find a sta- 
ble boson star configuration whose accretion disk shows 
the equivalent spectrum for a disk around a black hole of 
the same mass. 

We have also shown that for each value of the self- 
interaction coefficient A, once the mimicker boson star 
has been found for a given astrophysical combination 
with mass of the BHC and accretion rate, that same 
boson star configuration is the mimicker for other as- 
trophysical parameters combination. 

Most importantly, we also showed the way a boson star 
could be distinguished from a black hole by studying the 
light deflection. Considering that the black hole mimicker 
property of boson stars can be extended to other accre- 
tion disk models, it is possible that the light deflection 
could be together with the gravitational wave emission 
[24], HU , one of the few experiments that could discard or 
confirm the existence of boson stars. 



VI. CONCLUSIONS 

We have shown that stable spherically symmetric bo- 
son stars are black hole mimickers when the emission 
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